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1. Introduction 


The measurement of the transverse momentum and rapidity distributions of final-state par¬ 
ticles produced in the hadron-hadron and nucleus-nucleus collisions is considered to provide im¬ 
portant information about the collision process and the particle production mechanism [|^. The 
transverse momentum distributions of identified hadrons are fhe mosf common fools used fo inves- 
figafe sfrong inferacfions leading fo fhe creation of particles in fhe high-energy collisions due fo 
fhe momenfum balance in fhe plane perpendicular fo fhe beam direction. The missing fransverse 
energy is a key quantify in many searches for physics beyond fhe Sfandard Model, such as super- 
symmefry, exfra dimensions and dark maffer, and in fhe sfudy of fhe Higgs boson, neufrinos, W 
bosons and fop quarks 0. 

In fhe modern paradigm of science all known elemenfary particles are governed by fhe Bose- 
Einsfein or Fermi-Dirac sfafislics, i.e., fhey are bosons or fermions. The noninferacfing bosons 
and fermions given in fhe framework of fhe Bolfzmann-Gibbs sfafislics under fhe conslrainls of 
fhe grand canonical ensemble are described by fhe Bose-Einslein and Fermi-Dirac dislribulions, 
/ = l/[exp((e —/r)/r) ± 1]. These laws are followed by fhe photons irradiated in atomic and 
nuclear phenomena and by fhe parficles emitted in fhe hadron-hadron and heavy-ion collisions al 
low and infermediafe energies [i- However, fhe fransverse momenfum dislribulions of fhe parficles 
produced in fhe pp and AA collisions al EHC and RHIC energies can nol be described by fhe Bose- 
Einsfein and fhe Fermi-Dirac dislribulions of fhe Bolfzmann-Gibbs sfafislics given in fhe grand 
canonical ensemble. Thus, if we continue to accepl lhal fhe parficles af high-energies are bosons or 
fermions described by fhe Bose-Einslein or fhe Fermi-Dirac sfafislics, respeclively, Ihen we should 
find causes of changes of Iheir dislribulion funclions. 

The main purpose of fhis paper is fo verify fhe sfafislics of identified hadrons produced in 
pp collisions al high energies and to idenfify possible scenarios for fhe creation of hard hadrons 
wilh a large fransverse momenfum. Thereby, we invesfigafe fhe dala using fhe Bose-Einslein and 
Fermi-Dirac dislribulions of fhe Tsallis-facforized sfafislics and inlroducing fhe Iwo-componenl 
Bose-Einslein and Fermi-Dirac dislribulions of fhe ordinary Bolfzmann-Gibbs sfafislics. 

In fhe analysis of fhe EHC and RHIC dala, a Tsallis-like dislribulion gives salisfaclory fils fo 
fhe fransverse momenfum dislribulions However, hilherto in mosf works fhe Irans- 

verse momenfum dislribulions were fitted only by fhe classical Maxwell-Bollzmann dislribulion in 
both the Tsallis-factorized statistics ||^, 1^ 1^ 1^, 1^ and the two sources Boltzmann-Gibbs 

statistics |]J, 0 - 

The structure of the paper is as follows. In Section Q we briefly define fhe main formulas 
of fhe Tsallis sfafislics in fhe grand canonical ensemble. In Section ^ we remind fhe fransverse 
momenfum and rapidity dislribulions in fhe Tsallis-facforized sfafislics. The analysis of fhe dala 
and fhe resulls are discussed in Seclion The main conclusions are summarized in fhe final secfion. 


2. Tsallis statistics in grand canonical ensemble 

Eel us consider fhe equilibrium slalislical ensemble of fhe dynamical systems af fhe conslanl 
lemperalure T, volume V, chemical potential p and fhe Ihermodynamic coordinate z in a Ihermal 
conlacl wilh a heal balh. The system exchanges fhe energy and fhe matter wilh ifs surroundings. In 


2 




Non-extensive statistics effects in transverse momentum spectra of hadrons 


A.S. Parvan 


order to determine the probabilities of microstates, we consider the Tsallis statistical entropy [ p^ 
which is a function of the parameter q and a functional of the probing probabilities pp. 

S = (2-1) 

i ^ ^ i 

where kg is the Boltzmann constant, z = l/(^ — 1) and ^ G R is a real parameter taking values 

0 < q < °°. In the limit ^ > 1 the entropy (^) recovers the Boltzmann-Gibbs entropy. So = 

—ksY^Pilnpi. The probabilities of microstates are normalized to unity: 
i 

£p, = l. (2.2) 

i 

The parameter z takes the values —oo <z<—lforO<^<l and 0 < z < °° for 1 < ^ < oo. The 
expectation values of the Hamiltonian and the number of particles can be written as 


{H) = Y^piEi, 

i 

{N) = l^piNi, 


(2.3) 

(2.4) 


where Ei and Ni are the energy and number of particles in the /-th microstate of the system. In the 
state of thermal equilibrium the probabilities of microstates of the grand canonical ensemble can 
be found from the maximum entropy principle by the fundamental equation of thermodynamics at 
the fixed values of the variables of state {T,V,z,I-L) or by the Jaynes principle: 


Pi = 


1 + 


1=1 


1 A — Ei -\- pMi 
z+1 ksT 

1 A-Ei + pNP^ 


1 + 


z+1 


kBT 


(2.5) 

(2.6) 


where A is determined from Eq. (5^ and it is a function of the state variables, A = A{T,V,z,N). 
In the limit ^ 1 (z —)• ±°°) the probability p\ = exp[(flG — Ei + ixNi)/kBT] and the function A 

recovers the grand potential Q-g = —kBElnZc, where Zg = l^,exp[—(E, —/rA^,)/A:Br]. The expec¬ 
tation value (A) of the operator A can be defined as follows: 


(a) = La,- 

i 


^ 1 A-Ej + pNj 

^z+1 kBT 


(2.7) 


In fhe Tsallis sfafisfics, fhe mean occupafion numbers for fhe Bose-Einsfen, Eermi-Dirac and 
Maxwell-Bolfzmann sfafisfics of noninferacfing parficles in fhe grand canonical ensemble can be 
defined as 


{n 




= L G{npc 

{^pa} 


}n 


pG 


1 + 


1 A J^patipai^p P) 


Z+1 


kBT 


( 2 . 8 ) 


where ep = s/p^ + tfP- is fhe particle energy, G{npc,} = 1 for fhe Eermi-Dirac sfafisfics {n 
and fhe Bose-Einsfein sfafisfics (ns^ = 0) 1) • • •)°°) of parficles and G{npff} = 1/0 


pa 
pa '^pa • 


Maxwell-Bolfzmann sfafisfics {np(y = 0,1,... ,oo) of parficles |18]. 


= 0 , 1 ) 

for fhe 
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The transverse momentum and rapidity distributions of particles can be defined as 

In 


d^N V 


dprdy h? 


= 73 / d(ppT£pY,{t 


po! 


(2.9) 


0 


and 


^_V_ 

dy hr" 


In 


Jdq)J dpTPT£pY,{ 


npa), (2.10) 


0 0 

where pj and y are the transverse momentum and rapidity, respectively, and {np^j) are the mean 


occupation numbers given by Eq. (2.8). 


3. Tsallis-factorized statistics 


The Tsallis-factorized statistics is defined by the generalized entropy for the ideal gas in the 


grand canonical ensemble []13[]. In the Boltzmann-Gibbs statistics the entropy of the ideal gas in 
the grand canonical ensemble can be written as 


% — 

p 


{np)Gln{np)G + ;^(1 - T](n^)G)ln(l - T] {np)G) 


(3.1) 


where {np)G = l/(exp[j 8 (e^ “M)] + 0)> 0 = “1 for tho Bose-Einstein statistics, T] = 0 for the 
Maxwell-Boltzmann statistics and 77 = 1 for the Eermi-Dirac statistics of particles. In the Tsallis- 
factorized statistics, the Boltzmann-Gibbs entropy ( |3.1[ ) for the Bose-Einstein and the Eermi-Dirac 
statistics is rewritten, using the Tsallis prescription, in the following form [|T^, ^,'H] 


s = -gl, 

p 


{npY\nq{np) -f ;^(1 - r]{np)f\nq{l -r]{np)) 


(3.2) 


where 


In^(x) 


\-q 


(3.3) 


and <7 is a real parameter, q > Q. 
entropy (^T|) is rewritten as [ 22 ] 


Eor the Maxwell-Boltzmann statistics, the Boltzmann-Gibbs 


S = [{np)‘i\nq{np) - {np)] . 

p 


(3.4) 


The expectation value of the Hamiltonian and the mean number of particles of the system in the 
Tsallis-factorized statistics are defined as [13] 


{H) 

= sY^£p{npY^ 

p 

(3.5) 

(N) 


(3.6) 


p 
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where Ei and Ni are the energy and the number of particles in the /-th microstate of the system. 


Maximization of the entropies ( |3.2|) and ( |3.4i ) under the constraints (|3.5|) and ( |3.6|) leads to the 


mean occupation numbers for the Tsallis-factorized statistics []13|] 


{np) = 


1 


[l + {q-l)p{ep-n)]^ + T] 


(3.7) 


where T] = — 1 for the Bose-Einstein statistics, 77 = 0 for the Maxwell-Boltzmann statistics and 
T] = 1 for the Fermi-Dirac statistics of particles. Then the transverse momentum and rapidity 


distributions of particles can be written as [13] 


and 


d^N 


_ ^ ^ 

dpjdy J 
0 


271 


d(ppT£p{npY 


pjnijco^hy 


(271)2 ^[i_p(^_i)^(^ycoshy-/i)]‘?-i +riy 


(3.8) 


dN 

dy 


h? 


iTt 


dtp JdpTPT£p{npY 


gV 


PrdpT 


nij coshy 


(^[l + (^-l)j3(mr coshy-;U)]‘?-i +T]^ 


<1 ’ 


(3.9) 


where mj = YPt + Note that throughout the paper we use the system of natural units, h = 
c = kB = l. 


4. Analysis of transverse momentum distributions 


Let us verify that bosons and fermions produced in the pp collisions at LHC energies are 
described by the corresponding Bose-Einstein and Fermi-Dirac functions of both the Boltzmann- 
Gibbs statistics and the Tsallis-factorized statistics. We should check the fact that the particles 
created in the violent high-energy collisions are indeed bosons or fermions as those produced at 
lower energies. 

The transverse momentum distributions for the bosons (rj = —1), fermions (rj = 1) and Boltz¬ 
mann particles (rj = 0) in the Boltzmann-Gibbs statistics are given by 


d^N 

dpjdy 


gV 

(27r)2 


PT-mr coshy 


1 

nij coshv—^ 

e r +T] 


(4.1) 


These distributions correspond to the relativistic Bose-Einstein, Fermi-Dirac or Maxwell-Boltzmann 
ideal gases composed of many randomly moving identical point particles in equilibrium, which do 
not interact and which are given in the grand canonical ensemble at temperature T and chemical po¬ 
tential p in the volume V with imaginary walls. The chemical potential of antiparticles in Eq. (p~T|) 
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is taken with the opposite sign. At high-energies of the LHC the baryon chemical potential is zero, 
i.e. n =0. 

The creation of particles (bosons and fermions) in the high-energy pp collisions is a dynami¬ 
cal process. However, the transverse momentum distributions ( |2.9| ), (3.8) and (|4.l| ) of the Tsallis, 
Tsallis-factorized and the Boltzmann-Gibbs statistics can be applied only for the equilibrium sys¬ 
tems because they were obtained from the maximum entropy principle. Thus, the description of 
particle production process in high-energy collisions by these distributions leads to the substitution 
of the complex dynamical system with an equilibrium system of ideal gases. Is it possible to de¬ 
scribe such a complex dynamic process of creation of many particles only with a one equilibrium 
state of the system? Do all the particles produced in high-energy collisions originate from the same 
macroscopic state of the system or from different macroscopic states corresponding to different 
intervals of time of collision? 

Let us introduce the two-component transverse momentum distributions for the bosons, fermions 
and boltzmann particles given in the framework of the Boltzmann-Gibbs statistics 

d^N 


_ gVi 

dprdy (27r)" 


pjniT coshy 


1 

mj coshy— 

e +ri 


gV2 

+ coshy 


1 

nij’ coshy—fl2 

e ^2 -I- T] 


(4.2) 


This distribution corresponds to the ideal gas in two different macrostates of the equilibrium system 
(two sources). The first term describes the particles created in the macrostate {Ti,V\,pL\) and the 
second term is related to the production of the same type of particles from the other macrostate 
(22; V 2 ,jLt 2 )- The parameters {T,V,p) evolve during the evolution of the system. 

Let us analyze the transverse momentum spectra of particles created in the pp collisions at 
LHC energies. Figure 1 represents the transverse momentum distributions of pions kaons K^, 
protons and antiprotons produced in the pp collisions at energy = 0.9 TeV. The symbols rep¬ 
resent the experimental data of the ALICE Collaboration |0]. The solid and the dashed curves 
are the fits of the experimental data to the Tsallis-factorized function (^) and the two-component 
Boltzmann-Gibbs function (p~^, respectively, at rapidity y = 0. The values of the variables of state 
for the macrostates of the system for the Tsallis-factorized statistics and for the Boltzmann-Gibbs 
statistics are given in Tables 1 and 2, respectively. The pions and kaons follow the Bose-Einstein 
statistics. The protons and antiprotons follow the Eermi-Dirac statistics. The Bose-Einstein and 
Eermi-Dirac transverse momentum functions of the Tsallis-factorized statistics describe the exper¬ 
imental data very well. The values of the Tsallis-factorized temperature T, radius R of the system 
and the parameter q fluctuate in the dependence on the type of a hadron. The proton tempera¬ 
ture is much lower than the temperature of other hadrons; however, the volume of the system for 
protons is much larger than the volume of the system for other hadrons. Nevertheless, the value 
of the parameter q for protons is consistent with the values of the parameter q of other types of 
hadrons. In the pp collisions, the Tsallis-factorized volume of the system for all types of hadrons 
considerably overestimates the realistic volume of the system of two protons. The two-component 
Bose-Einstein and Eermi-Dirac transverse momentum functions of the Boltzmann-Gibbs statistics 
describe very well only kaons, protons and antiprotons. But the high-pr pions failed to be described 
by these two-component distributions. Eor all types of hadrons the two-component distribution is 
defined by two different sets of parameters {T,R) corresponding to two macroscopic states of the 
dynamic system. The first set of parameters corresponding to the lower part of pj--spectra (soft) 
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Pt (GeV) (GeV) 

Figure 1: (Color online) Transverse momentum distributions of identified hadrons produced in the pp 
collisions as obtained by the ALICE Collaboration 0] at ffs — 0.9 TeV. The solid and dashed curves are the 
fits of the data to the Tsallis-factorized distributions and the two-component Boltzmann-Gibbs distributions, 
respectively, at rapidity y = 0. 


is characterized by smaller values of the temperature and larger values of the volume in compari¬ 
son with the same parameters from the second set. The second set of parameters corresponding to 
the higher part of the pr-spectra (hard) is characterized by high temperatures and small volumes. 
The radius R\ is consistent with the radius of the system of two protons; however, the radius R 2 
corresponds to a strongly compressed system of two protons (core). Thus, in the pp collisions the 
high-P t hadrons are created in a system of a very small volume and large temperature. However, 
the low-pr hadrons are produced in a system of larger volume and smaller temperature. The Gibbs 
temperature T\ is about twice the Tsallis-factorized temperature T and the temperature T 2 is about 
twice the temperature T\. 


Figure 2 represents the transverse momentum distributions of the identified hadrons produced 
in the pp collisions at energy ^ = 0.9 TeV. The symbols represent the experimental data of the 
ALICE Collaboration |^]. The solid curves are the fits of the experimental data to the Tsallis- 
factorized distributions (^) and the dashed curves are the fits of the experimental data to the 
two-component Bose-Einstein and Eermi-Dirac transverse momentum functions of the Boltzmann- 
Gibbs statistics(p~^ at rapidity y = 0. The values of the variables of state for the Tsallis-factorized 
statistics and for the Boltzmann-Gibbs statistics are given in Tables 1 and 2, respectively. The 
mesons and (j) follow the Bose-Einstein statistics and the strange baryons follow the Eermi-Dirac 
statistics. The Bose-Einstein and Eermi-Dirac transverse momentum distributions of the Tsallis- 
factorized statistics describe the experimental data very well. The Tsallis-factorized temperature T 
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Px (GeV) 

Figure 2: (Color online) Transverse momentum distributions of identified hadrons produced in the pp 
collisions as obtained by the ALICE Collaboration at y/s = 0.9 TeV. The solid and dashed curves are the 
fits of the data to the Tsallis-factorized distributions and the two-component Boltzmann-Gibbs distributions, 
respectively, at rapidity y = 0. 


and radius R of the system strongly fluctuate in the dependence on the hadron type. The value of the 
parameter q fluctuates more slightly. The temperature of the strange baryon A differs essentially 
from the temperature of its antiparticle. The Tsallis-factorized volume of the system for K^, A and 
+ E^ considerably disagrees with the realistic volume of the system of two protons. The two- 
component Bose-Einstein and Fermi-Dirac transverse momentum distributions of the Boltzmann- 
Gibbs statistics describe the experimental data for all these identified hadrons very well. For the 
Boltzmann-Gibbs distributions the values of the temperature and the radius of the system also 
fluctuate considerably in the dependence on the type of a hadron. The first temperature Ti is smaller 
than the second temperature T 2 and the first radius R\ is larger than the second radius R 2 . The first 
set of parameters given by the small value of temperature and the large value of volume describes 
the low values of the transverse momentum spectra; however, the second set of parameters given 
by the large value of temperature and the very small value of volume describes the high values 
of the transverse momentum spectra. The radius R\ is consistent with the radius of the system of 
two protons and the radius R 2 corresponds to a compressed system of two protons. From the two- 
component Boltzmann-Gibbs distributions it is seen that the high-pr hadrons are produced from a 
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Type 

T, MeV 

R, fm 

q 

zVv 

Comments 

71+ 

75.94±0.57 

4.397±0.032 

1.142±0.001 

0.49 

Fig. 1 

71^ 

78.34±0.52 

4.262±0.027 

1.137±0.001 

0.39 


K+ 

58.32±6.22 

4.413±0.608 

1.172±0.007 

0.55 


K 

62.75±4.27 

4.070±0.348 

1.164±0.005 

0.30 


P 

33.44±13.2 

16.365±12.103 

1.147±0.009 

0.66 


P 

62.68±13.3 

5.449±1.991 

1.123±0.009 

0.65 



69.86±6.38 

3.531±0.441 

1.158±0.006 

0.78 

Fig. 2 

A 

104.79±34.00 

2.367±1.379 

1.096±0.017 

1.59 


A 

74.99±24.5 

4.273±2.688 

1 .111±0.012 

0.62 


0 

96 

1.228±0.104 

1.136±0.011 

0.60 


s-±r 

60 

3.221±0.176 

1.13 

2.21 



Table 1: Values of the parameters corresponding to the curves of fitting data to the Tsallis-factorized func¬ 
tions in Figs. 1-2. 


very compressed and heated system. The Gibbs temperature T\ is larger than the Tsallis-factorized 
temperature T and the Gibbs temperature T 2 is larger than the temperature T\. 


In Tables 1 and 2 we have the values of the variables of state for the Tsallis-factorized statistics 
and the Boltzmann-Gibbs statistics, respectively, for each type of hadrons produced in the pp 
collisions at y/s = 0.9 TeV measured by the ALICE Collaboration 0, ^|. The values of the 
parameters are not the same for different types of hadrons. However, at the fixed energy of collision 
the values of the parameters should be independent of the hadron type. Thus, let us find fhe mean 
values of fhe paramefers for fhe Tsallis-factorized sfafisfics and the Boltzmann-Gibbs statistics and 
compare them. For the experimental data of the ALICE Collaboration [Q, ^], the mean value of 
the Tsallis-factorized temperature T = 70.65 ±4.27 MeV, the mean value of the Tsallis-factorized 
radius of the system /? = 4.871 ± 1.151 fm and the mean value of the Tsallis-factorized parameter 
q = 1.138 ±0.003. For the two-component Bose-Einstein and Fermi-Dirac distributions of the 
Boltzmann-Gibbs statistics we have two sets of parameters. The mean value of the first temperature 
T\ = 148.38 ±4.59 MeV and the mean value of the first radius of the system Ri = 1.799 ±0.091 
fm. The mean value of the second temperature T 2 = 332.25 ± 12.15 MeV and the mean value of 
the second radius of the system R 2 = 0.419 ± 0.016 fm. 


It is seen that in the pp collisions at y/s = 0.9 TeV, the Boltzmann-Gibbs statistics suggests 
two regimes of production of particles. In the first regime, the Boltzmann-Gibbs soft hadrons are 
produced at the temperature Ti « 150 MeV from a heated spatial region of the volume consistent 
with the geometrical volume of two protons. However, in the second regime, the Boltzmann-Gibbs 
hard hadrons with a high transverse momentum originate from a hot superdense spatial region 
of the radius R 2 ~ 0.4 fm with the temperature T 2 ~ 330 MeV. In contrast, the Tsallis-factorized 
statistics presupposes that the particles in the pp collisions at y/s = 0.9 TeV are created from a 
large nonphysical volume of the radius « 5 fm at the small temperature T « 70 MeV. 
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Type 

ri,MeV 

Ri, fm 

r 2 , MeV 

R 2 , fm 

X^/v 

Comments 

71+ 

102.26±1.96 

3.424±0.059 

250.28±6.26 

0.934±0.044 

5.79 

Fig. 1 

71^ 

104.56±1.84 

3.348±0.052 

252.44±6.32 

0.905±0.044 

4.81 


K+ 

142.30±4.93 

1.671 ±0.077 

345.81±17.4 

0.380±0.035 

0.60 


K 

141.51±4.11 

1.677±0.060 

332.11±14.0 

0.399±0.032 

0.46 


P 

133.99±11.2 

2.216±0.432 

293.77±16.4 

0.420±0.058 

0.73 


P 

142.60±13.3 

1.860±0.366 

291.06±20.0 

0.416±0.074 

0.71 



141.12±5.10 

1.687±0.099 

340.85±8.77 

0.393±0.021 

0.75 

Fig. 2 

A 

179.17±38.5 

1.225±0.598 

352.36±45.4 

0.267±0.100 

1.69 


A 

170.38±26.3 

1.385±0.523 

348.92±32.0 

0.265±0.070 

0.81 


0 

191.58 

0.597±0.061 

473.51±116 

0.104±0.045 

0.66 


s-+r 

182.7 

0.700±0.112 

373.66 

0.129±0.013 

1.90 



Table 2: Values of the parameters corresponding to the curves of fitting data to the two-component 
Boltzmann-Gibbs functions in Figs. 1-2. 


5. Discussion and conclusions 

In conclusion, the Bose-Einstein and Fermi-Dirac distributions of the Tsallis-factorized and 
the two-component Boltzmann-Gibbs statistics lead to a satisfactory description of the data on the 
transverse momentum for the identified hadrons created in the high energy pp collisions. We have 
extracted the parameters T,V and q of the Tsallis-factorized statistics and the two sets of param¬ 
eters T\,V\ and 72, V 2 of the two-component Boltzmann-Gibbs statistics for each type of hadrons 
measured by the ALICE Collaboration [^, These parameters for different types of hadrons 
have different values for both the Tsallis-factorized statistics and the two-component Boltzmann- 
Gibbs statistics. The mean values of these parameters for the Tsallis-factorized statistics and the 
two-component Boltzmann-Gibbs statistics were obtained. The results of the Boltzmann-Gibbs 
statistics for the pp collisions at ^/s = 0.9 TeV suggest that the soft hadrons with low transverse 
momentum are created from the spatial volume comparable with the geometrical volume of two 
protons at the temperature of 150 MeV and the hard hadrons with high transverse momentum are 
produced from the very small volume in comparison with the volume of two protons at very high 
temperature of 330 MeV. In contrast, the volume of the system obtained by the Tsallis-factorized 
statistics is unusually large in comparison with the volume of two protons and the temperature is 
much less than both temperatures of the Boltzmann-Gibbs statistics. 

Note that the high-pr tail of the charged pion distributions cannot be described by the two- 
component Boltzmann-Gibbs statistics. This discrepancy for pions in hadronic collisions can be 
understood as a superposition of different production sources. 

Acknowledgments: This work was supported in part by the joint research project of JINR 
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References 

[1] W. Florkowski, Phenomenology of ultra-relativistic heavy-ion collisions (World Scientific Publishing 
Company, 2010) 


10 








Non-extensive statistics effects in transverse momentum spectra of hadrons 


A.S. Parvan 


[2] V. Khachatryan et al. (CMS Collaboration), arXiv: 1411.051 Ivl [physics.ins-det]. 

[3] K. Yagi, T. Hatsuda, Y. Miake, Quark-gluon plasma. From big bang to little bang (Cambridge 
University Press, 2008). 

[4] K. Aamodt et al. (ALICE Collaboration), Eur. Phys. J. C 71 (2011) 1655. 

[5] K. Aamodt et al. (ALICE Collaboration), Phys. Lett. B 693 (2010) 53. 

[6] K. Aamodt et al. (ALICE Collaboration), Phys. Rev. D 82 (2010) 052001. 

[7] G. Aad et al. (ATLAS Collaboration), New J. Phys. 13 (2011) 053033. 

[8] V. Khachatryan et al. (CMS Collaboration), Phys. Rev. Lett. 105 (2010) 022002. 

[9] M. Rybczynski, Z. Wlodarczyk, Eur. Phys. J. C 74 (2014) 2785. 

[10] J. Cleymans, G.I. Lykasov, A.S. Parvan, A.S. Sorin, O.V. Teryaev, D. Worku, Phys. Lett. B 723 (2013) 
351. 

[11] M.D. Azmi, J. Cleymans, J. Phys. G: Nucl. Part. Phys. 41 (2014) 065001. 

[12] J. Cleymans, D. Worku, J. Phys. G 39 (2012) 025006. 

[13] J. Cleymans, D. Worku, Eur. Phys. J. A 48 (2012) 160. 

[14] L. Marques, E. Andrade-II, A. Deppman, Phys. Rev. D 87 (2013) 114022. 

[15] P.-H. Liu, Y.-Q. Gao, B.-C. Li, Eur. Phys. J. A. 50 (2014) 123. 

[16] B.-C. Li, Y.-Z. Wang, P.-H. Liu, X.-J. Wen, Y.-E. Dong, Phys. Rev. D 89 (2014) 054014. 

[17] C. Tsallis, J. Stat. Phys. 52 (1988) 479. 

[18] A.S. Parvan, Theor. Math. Phys. 140 (2004) 977. 

[19] P. Buyukkilic, D. Demirhan, Phys. Lett. A 181 (1993) 24. 

[20] P. Pennini, A. Plastino, A.R. Plastino, Phys. Lett. A 208 (1995) 309. 

[21] A.M. Teweldeberhan, A.R. Plastino, H.G. Miller, Phys. Lett. A 343 (2004) 71. 

[22] T.S. Biro, E. Molnar, Phys. Rev. C 85 (2011) 024905. 

[23] K. Aamodt et al. (ALICE Collaboration), Eur. Phys. J. C 71 (2011) 1594. 


11 




